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Abstract. As an application of the new realistic three-dimensional (3D) formalism reported re- 
cently for three-nucleon (3N) bound states, an attempt is made to study the effect of three-nucleon 
forces (3NFs) in triton binding energy in a non partial wave (PW) approach. The spin-isospin depen- 
dent 3N Faddeev integral equations with the inclusion of 3NFs, which are formulated as function of 
vector Jacobi momenta, specifically the magnitudes of the momenta and the angle between them, 
are solved with Bonn-B and Tucson-Melbourne NN and 3N forces in operator forms which can be 
incorporated in our 3D formalism. The comparison with numerical results in both, novel 3D and 
standard PW schemes, shows that non PW calculations avoid the very involved angular momen- 
tum algebra occurring for the permutations and transformations and it is more efficient and less 
cumbersome for considering the 3NF. 
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INTRODUCTION 

As already known, the non-relativistic calculations of few-nucleon bound and scattering 
states are not able to reproduce experimental results for relevant observables such as 
the triton binding, even when considering various available realistic nucleon-nucleon 
(NN) interactions. Deviations can be attributed to wrong off-shell behavior of such 
potentials, relativistic corrections and, probably more important, 3NFs effect. In order 
to incorporate the 3NF corrections in a spin-isospin dependent 3D approach [1] for the 
triton, we have recently formulated the corresponding Faddeev equations in terms of 
the vector Jacobi momenta, specifically the magnitudes of the momenta and the angle 
between them, as well as the spin-isospin quantum numbers. 

We have shown that, for the full solution of the 3N bound system, the Tucson- 
Melbourne two-pion exchange 3NF can be included in a very simple manner in com- 
parison to the PW representation. As indicated in Ref. [2], according to the number 
of spin-isospin states that one takes into account, the formalism for both 3 H and 3 He 
bound states leads to only strictly finite number of coupled three dimensional integral 
equations, which at most for fully charge dependent case leads to 24 coupled equations. 
In this communication we present our numerical results for triton binding energy, ob- 
tained by solving spin-isospin dependent three-dimensional Faddeev integral equations, 
when considering the Bonn-B and Tucson-Melbourne NN and 3N forces. We would like 
to mention that our next task, which is currently underway, is to incorporate relativistic 



effects in 3N bound state calculations using the same 3D approach. 



FADDEEV EQUATIONS IN 3D REPRESENTATION WITH 3NFS 



In this section we briefly review the formalism of three-dimensional Faddeev integral 
equations in the realistic 3D approach. By considering the 3NF, the 3N bound state is 
described by the Faddeev equations 



| y/) = G tP\ y) + ( 1 + G t) G vgl |¥) , 



(1) 



where Gq = {E — Hq) is the free 3N propagator, the operator t = v + vGot is the NN 

transition matrix, P = PyiPn +^13^23 is permutation operator, the quantity V 12 3 defines 
the 3NF and \¥) = (1 + i D )|i//) is the total wave function. The representation of Eq. (1) 
in momentum space and in a non-PW scheme needs the following states in the 3D basis: 
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As shown in Fig. (1) the states of the 3D basis involve two standard Jacobi momentum 
vectors p and q. A comparison of basis states in both 3D and PW schemes show that 
in a standard PW representation the angular dependence leads to two orbital angular 
momentum quantum numbers, i.e., l\2 and ly. 
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whereas in 3D representation the angular dependence explicitly appears in the Jacobi 
vector variables. So it is clear that in 3D formalism there is not any coupling between 
the orbital angular momenta and corresponding spin quantum numbers. Therefore the 
spin quantum number of two-nucleon subsystem ^ and the third nucleon 53 = i couple 
to the total spin S and its third component M5. For the isospin quantum numbers, a 
similar coupling scheme leads to the total isospin T with its third component Mr- 
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FIGURE 1. 3N basis states in both 3D and PW representation. 

We would like to add the remark that in order to be able to evaluate the tran- 
sition and the permutation operators we need the free 3N basis states |pqy) = 
I p q m Sl m S2 m J3 m tl m t2 m t3 ) . To this aim when we are changing the 3N basis states | a ) 
to the free 3N basis states | 7) we need to calculate the Clebsch-Gordan coefficients 

(y\a) =g ya = (m Sl m S2 m S3 \(si2 ^)SM S ) (m h m t2 m t3 \(tn \)TM T ). 



By considering the symmetry property of the 3NF, the anti-symmetry property of the 
total wave function and the definition of anti- symmetrized NN /-matrix, i.e. a (\A)a = 
{\t{\ — Pn)\), the representation of Eq. (1) in Eq. (2) leads to: 
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As shown in Ref. [3], the evaluation of 3NF matrix elements (p qa |Vy 2 3l l I / ) f° r 2^- 
exchange TM force avoids the cumbersome nature of the PW representation and leads 
to simple expressions which are more convenient for numerical calculations. 



NUMERICAL RESULTS FOR 3 H BINDING ENERGY 

In this section we present our numerical results for the triton binding energy, obtained 
by solving the three-dimensional Faddeev integral equations (4). These coupled integral 
equations have been solved before in Ref. [2] without 3NF term. By solving eight 
coupled Faddeev equations for (j — j) spin-isospin states and by using the operator form 
of Bonn-B NN potential [4], our calculation for the triton binding energy converges to 
a value of E t = —8.152 MeV, whereas the PW calculations converges to E t = —8.14 
MeV for j™ x = 4. In order to be able to compare our numerical results with available 
PW calculations, we have used TM force [5] (with cutoff mass = 5.828 m K ) in 
an operator form [3] which is compatible with our 3D formalism. In table 1 we have 
shown the convergence of the triton binding energy in 3D approach as a function of the 
number of grid points. The corresponding PW results are also listed as function of total 
angular momentum j™*- ^ s demonstrated in this table, by using the Bonn-B and TM 
combination in the 3D approach the calculation of triton binding energy converges to 
E t = —9.75 MeV, whereas the corresponding result in a PW scheme for j^* = 1 yields 
E t = — 9.80 MeV. We should mention that this calculation is the first attempt toward the 
planned numerical investigations of 3 H binding energy with the most modern NN and 
3N forces, i.e., AV18 and modified Tucson-Melbourne three-nucleon force (TM'), and 
also consistent NN and 3N chiral forces. 



TABLE 1. The calculated triton binding energies E t of the three- 
dimensional Faddeev integral equations as function of the number of grid 
points in Jacobi momenta Nj ac and spherical angles N sp h, the number of 
grid points in polar angles is twenty. The corresponding PW results are 
listed for a comparison to our results. The used value for cutoff mass in 
TM 3NF is A n = 5.828/7%. 
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-8.15 [2] 
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